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I. INTRODUCTION 

Muon neutrino charged current quasielastic (CCQE) 
double differential cross sections have recently been mea- 
sured for the first time by the MiniBooNE collabora- 
tion The results demonstrate the inadequacy of the 
Relativistic Fermi Gas (RFG) as a model for the nuclear 
system. Indeed the RFG, presently used in the experi- 
mental analysis, underestimates the total cross section, 
unless an unusually large ad hoc value Ma = 1-35 GeV/c 2 
is used for the nucleon axial mass. 

The RFG model has the merit of treating properly the 
relativistic aspects of the problem. These cannot be ne- 
glected for the kinematics of MiniBooNE, where the neu- 
trino energy reaches values as high as 3 GeV. However, 
the RFG is clearly too crude to account for the nuclear 
dynamics, as is well known from comparisons with elec- 
tron scattering data. With this motivation several more 
sophisticated relativistic nuclear models have been ap- 
plied in recent years to neutrino reactions 043- The 
comparison with the new experimental data |l| allows 
for a detailed test of the corresponding predictions. 

Beyond the above-mentioned microscopic relativistic 
models, a phcnomcnological "SuperScaling" approach 
(indicated in what follows as "SuSA") has been pro- 
posed in (l5| . based on the assumed universality of the 
scaling function for electromagnetic and weak interac- 
tions. Analyses of inclusive (e, e') data have demon- 
strated that at energy transfers below the quasielastic 
(QE) peak superscaling is fulfilled rather well (l6l - fl8j 
(see also this means that the reduced cross sec- 

tion is largely independent of the momentum transfer 
(first-kind scaling) and nuclear target (second-kind scal- 



ing), when represented as a function of the appropriate 
scaling variable. From these analyses a phenomenological 
scaling function has been extracted from the longitudi- 
nal QE electron scattering response and used to predict 
neutrino-nucleus cross sections by assuming that this sin- 
gle universal scaling function is appropriate for all of the 
various responses involved (CC, CL, LL, T(VV), T(AA) 
and T'(VA); see [TH) and multiplying it by the corre- 
sponding elementary weak cross sections. 

Although being far more realistic than the RFG, the 
superscaling approach described above is based on some 
assumptions. First, in studies of inclusive QE electron 
scattering it assumes the equality of the longitudinal 
and transverse scaling functions. This property, which 
is known as scaling of the zeroth kind, has been tested 
in various models and shown to be violated to some ex- 
tent: for example Relativistic Mean Field (RMF) theory 
yields a transverse scaling function which is typically 20% 
or so larger than the longitudinal one [lU, [HI- In fact, 
this is exactly what is observed when one examines the 
existing L/T separated data. Once the effects that are 
expected to break scaling of the zeroth kind are removed, 
namely, inelastic contributions and effects stemming from 
two-particle-emission meson-exchange currents (see be- 
low) which are predominantly transverse in nature, one 
finds that the remaining transverse scaling function is 
clearly larger than the longitudinal one. Thus, when pro- 
ceeding to studies of CCQE cross sections, one should 
also expect to have some violations of scaling of the ze- 
roth kind as well. Second, the charged-current neutrino 
responses are purely isovector, whereas the electromag- 
netic ones contain both isoscalar and isovector compo- 
nents and the former involve axial-vector as well as vec- 
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tor responses. One then has to invoke a further kind of 
scaling, namely the independence of the scaling function 
of the choice of isospin channel — so-called scaling of 
the third kind. The interplay of scaling from the various 
contributions was first explored in [l3j . Finally, and most 
important, at energies above the QE peak scaling is vio- 
lated in the transverse channel by effects which go beyond 
the impulse approximation: inelastic scattering, meson- 
exchange currents (MEC) and the associated correlations 
which must be considered together with the MEC in or- 
der to conserve the electromagnetic current. 

In this paper we evaluate the double differential CCQE 
cross sections in the SuSA approach and discuss the im- 
pact of meson-exchange currents on the process. We are 
motivated by the fact that modeling at the level of the 
impulse approximation (as is the case for the RFG or for 
the spectral function approach of [20|-[22j]) under- predicts 
the measured CCQE cross sections and seems to call for 
a significant modification of the axial mass. However, 
more sophisticated approaches than the RFG such as 
SuSA and the other modeling discussed below are avail- 
able and the situation may not be quite so simple. For 
instance, previous non-relativistic calculations [23ll2~4| in- 
dicate that the 2p-2h excitations may be able to account 
for the large measured CCQE cross section, although a 
comparison of 2p-2h contributions with the MiniBooNE 
data is not yet available. However, it should be empha- 
sized that the kinematical regions explored under the 
integral over the neutrino flux extends over relativistic 
domains, and a relativistic treatment of the nuclear ex- 
citations is needed; this has clearly been shown to be 
necessary for electron scattering. In contrast, the 2p- 
2h MEC considered in the present work are taken from 
the fully relativistic model of [25J], where it was shown 
that relativistic effects are important to describe the nu- 
clear transverse response function for momentum trans- 
fers above 500 MeV/c. 



II. DOUBLE DIFFERENTIAL CCQE CROSS 
SECTIONS IN THE SUSA MODEL 



Following [15J we write the double differential CCQE 
neutrino cross section with respect to the muon scattering 
angle 9 and kinetic energy T M , at fixed neutrino energy 
E„, as 



dTf,d cos 9 



(Gf cos 9 c ) 2 k il 



£ u - — 



(1) 

where Gp is the Fermi constant, 9 C the Cabibbo an- 
gle, kfj, the muon momentum, Q 2 = uj 2 — q 2 the four- 
momentum transfer, with uj = E v — E^ and q = k v — k^, 
and 
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the nuclear response. The kinematical factors 1/j are 
given in (l5| : the nuclear response functions can be cast 
as 



qk F 
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where is the nucleon mass, kp is the Fermi momen- 
tum, Rf n - are the single nucleon responses, ip(q,u) is the 
RFG scaling variable (see, e.g., (2(| for its definition) and 
f(i/j) is the so-called superscaling function, containing the 
dependence on the nuclear model. In the RFG model the 
latter is a parabola limited to the region —l<ip<l: 
f($) = |(1 - ip 2 )9(l - V 2 )- In the SuSA approach, as 
already mentioned, it is given by a fit to the ex per imental 
longitudinal (e, e') reduced response function |19j . 

In order to compare with the MiniBooNE data we aver- 
age the cross section ((T|) over the neutrino energy flux pj . 
We use the Hoehler parametrization of the electromag- 
netic form factors and the value Ma = 1-03 GeV/c 2 
for the nucleon axial mass. In Fig. 1 we display the 
flux-integrated double differential cross section obtained 
within the SuSA model at fixed 9 as a function of the 
muon kinetic energy and compare with the MiniBooNE 
data [l[ . Since the experimental data are given in bins of 
cos 9, we display the results averaged over each angular 
bin. 

For most of the angle bins one sees that the SuSA 
results fall below the data and only for low scattering 
angles (for instance, look at the lowest two angle bins, 
0.9< cos 9 <1 and 0.8< cos6> <0.9) is the agreement rea- 
sonably good (note that an overall normalization error 
SN=10.7% should also be taken into account). However, 
one should be very cautious in applying models that are 
devised to work for quasi-free scattering, specifically the 
RFG model or the present SuSA approach. Such mod- 
els are not well suited to explaining the low-lying exci- 
tations in nuclei which arise from discrete states, giant 
resonances, etc. Indeed, when effects that fall under the 
heading "Pauli blocking" are tested in electron scatter- 
ing for excitations near threshold one does not see good 
agreement between experiment and such simple model- 
ing. The RFG model, the SuSA approach or any other 
models that lack the ability to address the complexity 
of the many-body problem in the near-threshold region 
are not supposed to be applied in this low q-uj regime. 
A proper treatment (for instance, using RPA with real- 
istic nuclear wave functions) of collective excitations is 
clearly required. Accordingly it is important to analyze 
how much of the integrated strength in the various panels 
in Fig. 1 arises from the low excitation region and how 
much from larger energies where the modeling may be 
more robust. 

To make such an assessment, in Fig. 2 we show the 
SuSA results obtained by integrating over the full neu- 
trino flux (solid lines, red online) and by artificially cut- 
ting the integral at energy transfer w=50 McV (dashed 
lines, green online). It clearly appears that at the most 
forward angles (upper left panel) when the results are cut 
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FIG. 1: (color online) Flux-integrated double differential cross section per target nucleon for the CCQE process on 12 C 
avaluated in the SuSA model and displayed versus the muon kinetic energy T M for various bins of cos 9. The data are from 
MiniBooNE The uncertainties do not include the overall normalization error <5iV=10.7%. 



at 50 MeV the cross section drops by about a factor of 
2, showing that roughly 1/2 of the cross section for such 
kinematics arises from the first 50 MeV of excitation. For 
more backward angles (upper right and lower panels) the 
cut effect is much weaker, indicating that low excitation 
energies arc only significant for the first angle bin. We are 
thus forced to conclude that approaches such as SuSA in 
the present paper (or the RFG, even when Pauli blocking 
effects are incorporated) should not be trusted for these 
very forward, significantly low-energy kinematics. 



In passing one should note that an additional correc- 
tion arises from the distortion of the outgoing muon wave 
function in the Coulomb field of the recoiling nucleus. A 
rigorous description of this effect is somewhat compli- 
cated as for the outgoing lepton it requires the use 
of distorted waves which arc eigenfunctions of the nuclear 
Coulomb field; however its main effects can be described 
using approximate approaches. In particular, the effec- 
tive momentum approximation 0, uBl |28| has been suc- 
cessfully applied to the case of medium-to-high energy 
leptons and is employed here. For muon kinetic energies 
T M > 200 MeV, as considered in the MiniBooNE exper- 
iment, Coulomb effects are below ~ 2%. Obviously, for 
smaller X^-valucs Coulomb distortion produces bigger ef- 
fects, these being on the order of ~ 10% for 7), ~ 90 — 100 
MeV (see @]). 
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FIG. 2: (color online) solid lines (red online): flux-integrated 
cross sections calculated in the SuSA model for specific bins 
bin of scattering angles, dashed lines (green online): a lower 
cut u — 50 MeV is set in the integral over the neutrino flux. 



III. TWO-PARTICLE TWO-HOLE 
MESON-EXCHANGE CURRENTS 

In this section we evaluate the contribution of meson- 
exchange currents to the CCQE cross section. The MEC 
are two-body currents, and therefore can excite both one- 
particle one- hole (lp-lh) and two- particle two-hole (2p- 
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FIG. 3: (color online) Same as Fig. 1 for the SuSA model, but now including 2p2h meson-exchange currents. 



2h) final states. 

Most of MEC studies of electromagnetic (e,e') pro- 
cesses performed for low-to-intermediate momentum 
transfers in the lp-lh sector (see, e.g., [2914321 ]) have 
shown a small reduction of the total response at the 
quasielastic peak, mainly due to diagrams involving the 
electroexcitation of the A resonance. These roughly com- 
pensate the positive contribution due to correlation dia- 
grams, where the virtual photon couples to a correlated 
pair of nuclcons. In the present work we shall therefore 
neglect them and restrict our attention to 2p-2h final 
states. 

The impact of pionic 2p-2h MEC on inclusive electron 
scattering reactions has first been evaluated in the RFG 
framework in [33j |. where a non-relativistic reduction of 
the currents was performed. Fully relativistic calcula- 
tions have been developed more recently in [gjj HJ] . It 
has been found that the MEC give a significant positive 
contribution to the cross section, which helps to account 
for the discrepancy between theory and experiment in the 
"dip" region between the quasielastic and A-resonance 
region. Moreover, the MEC have been shown to break 
scaling of both first and second kinds [H[ . 

In this paper we use the fully relativistic model of 
[25I ] . where all many-body diagrams containing two pi- 
onic lines that contribute to the electromagnetic 2p-2h 
transverse response were taken into account. Note that 
in lowest order these affect only the transverse polar vec- 
tor response, R)£ v . As shown in Fig. 3, the 2p-2h MEC 
tend to increase the cross section, yielding reasonable 
agreement with the data out to cos 6 ~ 0.6. At larger 



angles the disagreement with the experiment becomes 
more and more significant, and the meson-exchange cur- 
rents are not sufficient to account for the discrepancy. 
We also note that recent studies using somewhat differ- 
ent assumptions [36| yield similar results for 2p-2h MEC 
contributions, although these tend to be somewhat larger 
than the earlier results employed in the present work, 
and would lead to somewhat better agreement with the 
CCQE data at larger angles. 

The same conclusion discussed above can be drawn 
by plotting the cross section versus the scattering angle 
at fixed (Fig. 4): the inclusion of two-body currents 
improves the agreement with the data at low scattering 
angles, but some strength is missing at higher angles, 
especially for low muon momenta. 

Again some caution should be expressed before draw- 
ing definitive conclusions from the agreements or dis- 
agreements seen in Figs. 3 and 4. For instance, as al- 
ready mentioned, there are strong indications from RMF 
studies as well as from QE (e, e') data that scaling of 
the zeroth kind is only approximate and that the vector 
transverse response should be enhanced over the strict 
SuSA strategy employed in the present work. Moreover, 
in different language, namely that of extended RFG mod- 
cling where lp-lh and 2p-2h excitations are incorporated, 
a fully consistent treatment should take into account not 
only the MEC contributions of the present study but also 
the correlation diagrams that are necessary in order to 
preserve the gauge invariance of the theory. In an in- 
finite system like the RFG these diagrams give rise to 
divergencies which need to be regularized. A treatment 
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FIG. 4: (color online) Double differential z/ M CCQE cross section for 12 C integrated over neutrino flux versus the outgoing muon 
scattering angle for various bins of the muon kinetic energy T M . Solid lines (red online): SuSA; dashed lines (green online): 
SuSA model including the MEC 2p-2h contribution. In the light of the discussion of Fig. 2 given above, results are given only 
for the region cosf <0.9. 



of these contributions has been performed recently for 
(e, e') in (36j . where they were shown to be of the same 
order as the MEC. The model of [36[ contains similar in- 
gredients in the treatment of 2p-2h excitations to those 
in [23|, Hil with the addition of being fully relativistic. 
These contributions also enhance the cross sections be- 
yond the results shown above and so might be respon- 
sible for the residual disagreement. Indeed the trend is 
encouraging: their effect grows with increasing momen- 
tum transfer and hence they have a greater impact for 
large neutrino-muon angles than for the forward direc- 
tion where the results shown in Figs. 3 and 4 are already 
reasonably successful. It is too early to say for sure, how- 
ever, since there is at present no completely consistent 
relativistic model that is capable of incorporating all of 
the effects discussed above. 



IV. CONCLUSIONS 

In summary, we have applied the phcnomenological 
model based on electron scattering data, elaborated in 
[Til , to CCQE neutrino reactions and compared the 
results with the recent MiniBooNE double differential 
quasielastic cross section data at all available kinemat- 
ics. In addition to presenting detailed results for the 
so-called SuSA approach, two specific issues have been 
addressed in the present work: (1) the role played by 2p- 
2h MEC contributions has now been explored, and (2) 
cross sections at small angles have clearly been shown 
to be related to the regime where low-energy nuclear ex- 
citations dominate and thus where quasi-free modeling 
must be viewed with suspicion. The strict SuSA predic- 
tions show a systematic discrepancy between data and 
theory, where they tend to underestimate the data es- 
pecially at large muon scattering angles and low muon 



energies. When 2p-2h MEC contributions are included 
the situation is different: inclusion of the 2p-2h contribu- 
tions yields results that are compatible with the data for 
9 < 50° (excluding the most forward angles where quasi- 
free modeling must be questioned, as stated above), but 
lie below the data at larger angles where the predicted 
cross sections are smaller. These two-body currents arise 
from microscopic relativistic modeling performed for in- 
clusive electron scattering reactions and they are known 
to result in a significant increase in the vector-vector 
transverse response function, in concert with QE elec- 
tron scattering data. It should, however, be remembered 
that the present approach still lacks the contributions 
from the correlation diagrams associated with the MEC 
which are required by gauge invariancc; these might im- 
prove the agreement with the data, as suggested by the 
results of [36| for inclusive electron scattering. Finally, 
we note that alternative approaches such as relativistic 
mean field theory also lead us to expect an enhancement 
over the results shown in the present work, as, in fact, 
are observed for QE electron scattering data in a similar 
kinematic region. Work is in progress to resolve several 
of these issues, specifically, to perform a detailed study 
of modeling versus experiment for inclusive QE electron 
scattering, to extend the analysis based on RMF both 
for electron scattering and for neutrino reactions, and to 
incorporate the missing pieces mentioned above that are 
required to restore gauge invariance. Once these are in 
hand it will be appropriate to re- visit the issue of the 
anomalous axial mass. 
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